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I. 



INTRODUCTION 



For many practical applications one needs accurate val- 
ues of thermodynamic properties of solids and liquids 
both at low pressures, where direct measurements are 
available, and in regions of temperature and pressure 
where data are absent. This requires a model Hamil- 
tonian for the system which gives accurate results for the 
Helmholtz free energy. For a metal crystal the free energy 
consists of three contributions, 



$ {V) + F E (V,T) + F T (V,T) 
$ (V) + F e i(V,T) + F ep {V,T) 
+ F ph (V,T) + F anh (V,T). 



(1) 



§o(V), the static lattice potential, represents the total 
energy when the ions are located at lattice sites and the 
electrons are in their ground state. Fi(V,T) is the free 
energy from ion vibrations, and Fe(V,T) is the free en- 
ergy associated with thermal excitation of electrons from 
their ground state. Fj(V, T) is the dominant temperature 
contribution. It consists of the quasi-harmonic phonon 
contribution, F p h{V,T), plus the small anharmonic term 
F an h{V, T) which expresses phonon-phonon interactions. 
Fe(V,T) consists of F e i(V, T) representing the thermal 
excitation of independent electrons, plus F ep (V,T), the 
contribution from interactions between electronic excita- 
tions and phonons. 

In this paper we determine the magnitude and temper- 
ature dependence of the non-adiabatic contributions to 
the free energy from the electron-phonon interaction for 
different metals and compare with previous treatments. 

Our calculations are done for a constant density p to 
eliminate concern for the density dependence of phonon 
frequencies and electron-phonon interaction matrix el- 
ements. The density is that at the temperature T p , 



where the phonon frequencies arc measured. The melting 
temperature at this density is higher than the custom- 
ary zero-pressure melting temperature. Our calculations 
cover the range from T = to above T m . 

The paper is organized as follows: In section[n]we will 
introduce the expressions used for the electron-phonon 
part of the free energy. In section ITTTl our results are pre- 
sented with a more detailed discussion of the numerical 
methods used. 



II. ELECTRON-PHONON FREE ENERGY 

A. Analytic Form of the Free Energy 

Let us briefly discuss the physical origin of the con- 
tributions to F ep , and show in the process that our 
formulation is free of double-counting errors. In elec- 
tronic structure theory we calculate the electronic ground 
state energy with static nuclei, as a function of the nu- 
clear positions, {rx} = ?k=\,- - ,n- This energy is the 
ground state adiabatic potential <j>({rR-}), which is re- 
solved into $o + &ph + ®anh, and which together with 
the nuclear kinetic energy gives the free energy contri- 
bution $o + F p h + F an h in eq. Q. We next do the 
same calculation for the excited electronic states, labeled 
n = 1, 2, • • • of the energy E n ({rk-}) with static nuclei. 
Here only the excitation energies E n — $ arc considered 
since $ is already included in the nuclear motion con- 
tributions. Excited electronic states with nuclei at the 
crystal lattice sites are calculated in electronic structure 
theory, and their energy levels are expressed in the elec- 
tronic density of states, which provides the free energy 
contribution F e i in eq. JIJ. Next the thermally averaged 
vibrational contributions to the excited energies E n — $ 
are calculated, and these contributions yield the adia- 
batic part F ad of the electron-phonon free energy, eq. 



2 



© . Finally the non-adiabatic correction to all the elec- 
tronic energy levels, ground and excited states alike, is 
calculated by allowing the nuclear motion to mix the elec- 
tronic states. This produces the terms F" a + F^ a , eqs. 
and 10. The non-adiabatic part can formally be at- 
tributed to the terms in the trace of the partition function 
in which the nuclear kinetic energy operator operates on 
the electronic wavefunctions. A detailed discussion of 
the resolution of the crystal Hamiltonian, and the cor- 
responding free energy contributions may be found in 
Wallacei Sees. 4 and 18 and pp. 91 - 94. 

While this theory is valid for metals in general, an 
approximation is available from which we can calculate 
F ep from previously calibrated models for the nearly- 
free electron metals. This is pseudopotential perturba- 
tion theory, which starts from free electrons in zeroth or- 
der, and treats the screened electron-ion interaction (the 



pseudopotential) as a perturbation. Electronic band- 
structure effects then arise in standard perturbation the- 
ory. The electron-phonon theory becomes a double per- 
turbation expansion, in the pseudopotential and in the 
displacement of nuclei from equilibrium. But the dis- 
placement expansion converts to a pseudopotential ex- 
pansion, and the leading contributions to F ep are all of 
second order in the pseudopotential. The derivation may 
be found ini (Sec. 18). 

The electron-phonon contribution is written as a sum 
over the three terms mentioned in the previous para- 
graph. They are 

F ep = F ad + Fl ia + F 2 " Q , (2) 
where the single contributions are given by 
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Here and in the following we will calculate and quote 
our results as per atom, fp is the Fermi-Dirac distribu- 
tion function at finite temperature and gp is the same at 
T = 0. rit\ is the Bose-Einstein distribution function at 
finite temperature and fj^ x is the polarization vector of 
the phonon branch A for wave vector k which is inside the 
Brillouin zone. Q is a reciprocal lattice vector and is 
the frequency of a phonon mode. U(k + Q) is the Fourier 
transform of the pseudopotential for momentum transfer 
k+Q. In our calculations we use two models for the pseu- 
dopotential, the Harrison and the Ashcroft models, which 
are screened by exchange and electron-electron interac- 
tions (detailed forms are given in Wallace^, pp. 312). 
The pseudopotential parameters are listed in Table [H] 
ep is the electron energy measured relative to the Fermi 
energy. 

In the above formulas, every term is of second order in 
U(q). The appearance of free electron energies in the de- 
nominators results from the use of perturbation theory. 
The pseudopotential modification of the electron wave- 
functions gives rise to the last term in braces in eq. J3J). 
An alternate derivation, by means of Matsubara frequen- 



cies, is outlined in the appendix. 

This is the form of the electron-phonon contribution to 
the free energy originally calculated by Eliashberg— who 
studied the electron-phonon interaction in the supercon- 
ducting state. It is known that the non-adiabatic contri- 
bution dominates the adiabatic one at low temperaturai 
and we will focus on F™ m this paper. F ad will be ex- 
amined in a future publication. 



B. Analytic temperature dependences 

From eqs. Q and JSJ we estimate the low-T depen- 
dence of the two non-adiabatic contributions with the 
leading order terms of a Sommcrfeld expansion^ (pp. 45). 
This is certainly accurate for nearly free electron metals 
which have smooth density of states around the Fermi 
surface, but is less certain for transition metals and ac- 
tinides because their electronic density of states tends to 
fluctuate strongly around the Fermi surface which makes 
power series expansions around ep less accurate. 



In order to analyze the first non-adiabatic contribution 
we rewrite the Fermi-Dirac factor as follows: 



fp = 9p +(.fp- 9p) 



(0) 



At high temperatures we can expand the other tempera- 
ture dependent factor 



(3hw 
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720 



(7) 

in powers of T and get a linear T dependence in leading 
order. At low temperature, the factor eq. J7J) will go over 
to a constant and we find that the temperature depen- 
dence of the ground state part of F" a given by the first 
term in eq. © will be given by 



const, at T = 
kT at high T 



(8) 



The contribution from the term in parentheses in eq. (jSJ) 
is restricted in phase space to a small volume around the 
Fermi surface, whereas the ground state contribution is 
restricted only to p < pp. Since the (fp — gp) contribu- 
tion has much less phase space than the gp contribution, 
the ground state contribution will dominate at all tem- 
peratures and we will not sec any additional temperature 
dependence beyond eq. JjjJ. 

Considering the factor fku/ (^ep — epj in F™ a we expect 

this contribution to be of order K(u>) /ep smaller than 
F™ and we will be able to neglect it. We calculated F™ a 
for Na to verify these estimates and present it in section 



2. 



Taking a closer look at the integrand of the second 
non-adiabatic contribution we notice that all factors but 
one are positive for all points in the integration region. 
Only the energy difference in the denominator can have 
different signs depending on the momenta p and p' = 
p+k+Q. We estimate its sign by dividing the integration 
phase space into two regions 
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(Region 1) 
(Region 2) 



(9) 



At zero temperature due to the Fermi-Dirac factor, 
fp(l — fp^j, the width of the phase space for Region 
1 is twice the phonon energy, 2hjj, in energy around the 
Fermi surface and very small compared to Region 2 which 
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FIG. 1: w B -dependence of F?" for Na at T = 200 K. 
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FIG. 2: F? a for Na. 



extends everywhere outside Region 1. i^ a (T = 0) will 
therefore be positive. Increasing temperature will widen 
the phase space for Region 1 because of the softening of 
the Fermi surface and will increase its negative contri- 
bution to the total but will leave the contribution from 
Region 2 mostly unaffected. The overall effect will be to 
reduce the magnitude of F£ a , the more so with increasing 
temperature but we expect F™ to remain positive. 

From a Sommerfcld expansion of the Fermi-Dirac 
factor for small temperatures we estimate the low- 
temperature dependence of F™ to be given by 



jpna 
b 1 



C 2 + A 2 T 



(10) 



Previous work on the electron-phonon contribution to 
the free energy focused on calculating the specific heat 
and missed the factor C 2 which we will calculate in sec- 
tion IIIII An extensive survey of the theoretical calcula- 
tions is given by Grimvall^ (Tables III- VI). We will also 
calculate the factor A 2 and compare with experiment and 
previous estimates for different materials. 

As has been shown by previous authors, the low tem- 
perature specific heat of normal metals has a T 3 In T con- 
tribution which comes from electron electron interaction, 
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both those induced by phonons and those due to the 
Coulomb interaction. Coffey and Pethicki and Danino 
and Overhauses^ for instance derived this contribution 
for a Debye model. We are using an Einstein model in 
our calculations of the non-adiabatic parts and can not 
pick up the T 3 In T term since it depends on the existence 
of acoustic phonon modes. It is known however that this 
contribution is small compared to the others and we will 
not calculate it. 



III. 



RESULTS 



A. Einstein Approximation 

In this section we want to study the two non-adiabatic 
contributions using an Einstein model. We will try to 
answer two main questions: Will we get the correct tem- 
perature dependence from such a simplified model? How 
accurate will our results be? 

The details of the phonon model enter the single con- 
tributions through the factor 



k 



(11) 



in the integrands and through their explicit dependence 
on the phonon frequencies, w^ A . The former contains 
details of the phonon spectrum in terms of the polar- 
ization vectors of the phonon branches and turns out to 
have a fairly weak effect on the integral, whereas the ef- 
fect of the latter can be very strong for acoustic phonon 
branches. When k goes to zero at the zone center— the 
phonon frequencies of acoustic branches vanish. If the 
integrand plus volume factors of k 2 is well-behaved and 
finite for — > we can expect to be able to approxi- 
mate the integral using a suitably chosen Einstein mode. 
If the integrand plus volume factors diverges for vanish- 
ing wg A the details of the phonon spectrum will be im- 
portant and we can not expect to be able to approximate 
the full phonon spectrum with a single frequency. There 
might be a situation of course in which a single phonon 
frequency will result in the correct temperature depen- 
dence, but without calculating the integral using a real- 
istic phonon dispersion we have no way of knowing what 
that single frequency should be. 

Both non-adiabatic contributions, 
tor 



F™2, have the fac- 



1 



e P e p+k+C 



in common. At the zone center this turns into 

1 



(12) 



e P e p+k+Q 



(13) 
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FIG. 3: ^-dependence of F£ a for Na at T = 200 K. 
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FIG. 4: F£ a for Na. The inset shows the low temperature 
region. 



which is well behaved^. The additional factor 



lim hu %x ( n %x + | ) = k B T 



(14) 



in F™ a is also well behaved at the T point. We conclude 
that a single phonon frequency will reproduce F™?; with 
the correct temperature dependence. We also expect it 
to be accurate and we will quantify this last statement 
later in this paper. 

Figs. Hand El show the non-adiabatic contributions for 
an Einstein model at T = 200 K. Our a priori estimate 
for the Einstein frequency is lue = (w), where (lu) is 
the average over the Brillouin zone. For reference we 
included our estimates for the average phonon frequency 
in the graphs and in Table HJ As expected, the frequency 
dependence is very weak and we can neglect it. 



B. Numerical Techniques 

The two sums over the electron and phonon momenta, 
y~]p and yiggj arc over all of reciprocal space. Assuming 
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a large sample size we therefore have the option to con- 
vert one or both of these sums into integrals, using the 
well known relationship 



E 



NV A 
(2irf 



d A p, 



(15) 



where N is the number of ions in the sample, and Va the 
volume per ion. In general, due to the completeness of 
the phonon eigenvectors in 3 dimensional space, 



In an Einstein model the sum over phonon branches can 
be performed outside the integral since the only factors 
that depend on the branch index, A, arc contained in eq. 
(|16[) . This means that the integrand does not depend 
on the crystal structure. If we combine the double sum 
SfcQ with where p' = p + k + Q we can replace both 
sums, Y%-,, in the two non-adiabatic contributions with 
integrals without loss of generality. Eqs. Q and J5J) can 
then be written as 
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(18) 



where V stands for the principal part. We would like to 
mention that neither one of the two non-adiabatic contri- 
butions has any obvious upper limits for the integration 
on p' and we need to make sure that we understand the 
convergence behavior of their integrands for large p' . 

The pscudopotcntial, U(q), converges as q~ 2 for large q 
and we conclude that the second non-adiabatic contribu- 
tion, KT a , converges as (p') _4 . Due to the additional fac- 
tor of 1/ (e p — e p i) in F" a , the first non-adiabatic contri- 
bution converges as (p')~ 6 . We can therefore expect good 
convergence behavior for the two non-adiabatic contribu- 
tions in terms of the integration over p' = p + k + Q. Fig. 
[S] shows this nicely. p max is the upper integration limit 
oip'. 

The 6-dimensional integration of eqs. (|17fl and 118|) 
was done using the VEGAS methodSiiS. We took special 
care to ensure convergence around the divergent poles of 
the integrands when calculating the principal part. This 
was done by first reducing the poles from second to first 
order and then removing the divergence by subtracting 
a suitably chosen function from the integrand and in- 
tegrating its pole by hand. This effectively smoothed 
the integrand sufficiently so that our VEGAS calculation 
converged much more rapidly. It also improved the re- 
liability of the VEGAS error estimate and therefore the 
accuracy of our calculation. Since we are using spherical 
coordinates, this procedure has the advantage that we 
can very efficiently limit the integration region around 
the Fermi surface by using appropriate integration lim- 
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FIG. 5: Convergence of F?" and F 2 ™ for Na at T = 200 K. 



its. It is also simpler to implement than a tetrahedron 
method in which the Brillouin zone is divided into tetra- 
hedra which are then summed up. 



C. Non— Adiabatic Contribution to F ep for 
different materials 

We calculated F[ la for Na and F 2 " a for Na, K, Al, and 
Pb. Table|Ushows the material specific constants used in 
our calculations. 

Figs. [3 and 0] show our results for Na using a single 
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FIG. 6: Ff a for K, Harrison model. 

TABLE I: Material constants 

Na K Al Pb 
structure bcc bcc fee fee 
T p [K] 90 9 80 80 

p [g cm- 3 ] 1.005 0.904 2.731 11.55 
{hco) [meV] 10.53 6.42 25.79 5.90 
a [A] 4.23 5.24 4.05 4.95 

p F [A- 1 ] 0.92 0.74 1.75 1.57 
e F [eV] 3.24 2.11 11.66 9.45 
Tmeit [K] 407 368 1234 722 



phonon mode at Huje = 10.53 meV for F™- F\ a exhibits 
linear temperature dependence for higher temperatures 
and even at temperatures as low as 50 K we hardly sec de- 
viations from linear behavior. This confirms our previous 
estimate of this contribution's temperature dependence 
based on a phase space argument. In section Hi Bl we had 
argued that when we rewrite the Fermi-Dirac distribu- 
tion factor in eq. ijTTj) as eq. (JBJ we expect the (^factor 
to be dominant due to its much larger phase space. Since 
gp is not temperature dependent, the only temperature 
dependent factor is the phonon distribution which is lin- 
ear at high temperatures. As we had argued in the same 
section also, the magnitude of F™ a is much smaller than 
the magnitude of F£ a . Based on this magnitude differ- 



TABLE II: The pseudopotential parameters. Parameters for 
the Harrison model are p and /3, while the single parameter 
for the Ashcroft model is p. 

Na K Al Pb 



Reference - (p. 406) - (p. 406) - (p. 415) — (p. 502) 



z 


1 


1 


3 


4 


C 


1.81 


1.77 


1.90 


1.81 


p [a B ] 


0.50 


0.69 


0.24 




$ [% a%] 


37 


66 


47.5 




p [a B ] 






1.117 


1.12 
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q/2p F 

FIG. 7: Harrison and Ashcroft pseudopotential models for 
Al. The parameters are given in Table EH 



ence we will neglect F™ a from this point onward. 

Fig. 01 shows F-f a and as expected this contribution 
is positive all throughout the temperature regime shown. 
Starting from T = the curve has negative curvature but 
quickly switches to positive curvature and never crosses 
zero even at high temperatures. From the inset in Fig. 0] 
we get the values of Ci and A 2 listed in Table ITTT1 

Fig. shows our result for F™ for K. The magnitude 
of the F™ contribution is smaller by about a factor of 
3 compared to Na. The difference is presumably largely 
due to the difference in ep which affects the size of the 
pseudopotential U (q) . Since the ratio of Fermi energies is 
about 3/2 and the pseudopotential enters the free energy 
squared, we get about a factor of 2 difference between 
the two materials just based on this term. 

The pseudopotentials for Al are shown in Fig. Be- 
tween zero and q/(2pp) ~ 1 the two models are iden- 
tical. Between 1 < q/(2pp) < 4 the Harrison model is 
larger than the Ashcroft model and for qj (2pf) > 4 both 
models slowly approach zero. Fig. [S] shows our results 
for Al which we calculated for these two pseudopotential 
models. F% a using an Ashcroft pseudopotential model is 
shifted by about 0.05 meV to lower energies compared 
to the Harrison model. The qualitative temperature de- 
pendence is not affected. The difference between the two 
models in terms of their convergence to zero for large 
q therefore does not appear to be relevant in affecting 
the temperature dependence of F^ a . The difference in 
magnitude for intermediate q does enter the overall mag- 
nitude of F™ and makes a difference of about 20 - 30%. 
The inset Fig. [SJshows a fit using eq. ^| and the results 
for C'2 and A2 are listed in Table ITTT1 

Fig. |21 shows our results for Pb. In Fig. Q^wc show 
the entropy calculated from the free energy and compare 
our result with earlier results by Grimvall— , Fig. 5.19. 
Our results are scaled with an Einstein temperature of 
(hu) = 5.9 meV = 68.5 K. We find very good agree- 
ment with Grimvall's results. Grimvall calculated the 
total entropy, i.e. S GrimvaU = S ad + S[ la + S% a for Pb 
up to T ~ 1.4 Te- Up to this temperature, our purely 
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FIG. 8: F 2 a for Al for different pseudopotential models. 
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FIG. 9: F 2 ' M for Pb, Ashcroft pseudopotential model. 



non-adiabatic entropy, S™, agrees very nicely with his 
results^. The adiabatic contribution remains small com- 
pared to the non-adiabatic up to much higher tempera- 
tures. We shall address this in a future publication. 



IV. CONCLUSIONS 

A thorough analysis of the first non-adiabatic term, 
Fi a , reveals a ground state contribution, oc gp, and a 
contribution due to excited electronic states, oc {fp — gp), 
cf. eq. @). We find that the ground state term, since 
it's temperature dependence is determined by the phonon 
factor alone, leads to a T-dependcnce of F" a as given in 
eq. JHJ). Electronic excitations above the ground state, 
the (fp — gp) term, will lead to a temperature dependence 
which is determined by the leading terms in a Sommerfcld 
expansion and give a T 2 -dependence at low temperatures 
for Fi a . Using a phase space argument we found that 
the ground state term will dominate the temperature de- 
pendence and essentially determine it. The contribution 
from the excited electronic states to the temperature de- 
pendence will be much smaller than the ground state 
contribution and not visible in our calculations. We also 
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T[K] 

FIG. 10: FT and S 2 a {T)/ (-2 A 2 T), where A 2 was defined 
in eq. (I1UL compared with earlier results by Grimvall on Pb. 
F 2 a is plotted with axes on the left and bottom, S 2 a and 
GrimvalPs results are plotted with axes on the right and top. 



TABLE III: F 2 na (T = 0) and curvatures at low-T (Units: C 2 
[1(T 8 eV/atom], A 2 [1(U 5 eV/atom K 2 ]) 







Na 


K 


Al 


Pb 


C-2 


(Harrison) 


2.270 


0.8543 


32.09 




(Ashcroft) 






27.04 


6.839 




(expt. °) 


-0.10 


-0.13 


-0.19 


-0.86 


A 2 


£ (Table 27) 


-0.14 


-0.14 


-0.30 




(ppt ") 


-0.10 


-0.12 


-0.24 


-0.91 




(Harrison) 


-0.0878 


-0.0889 


-0.1316 






(Ashcroft) 






-0.1285 


-0.6039 



a A 2 is calculated from experiment using low-temperature specific 
heat data. From eq. we nn d C ep = -2A 2 T. At low T, 

Cexpt = r expt T and C th = T ba (1 + A) T. Our C ep = A r ba T and 
we find A 2 = — i (F exp t — Ti, s ). F exp t taken from Kittel™, Table 
on p. 157. The analysis was done by Allen™, Table I and eq. (6). 

b A 2 is calculated from pseudopotential perturbation theory, and 
given by A 2 = — ^r;, s A. A is taken from the data collection of 
Grimvall, Tables III, IV, and V. 



expect the magnitude of F" a to be oc {{fmj) /e_f) F 2 " a an d 
hence negligible compared to ■ Our numerical results 
confirm our suspicions, cf. Fig. [21 The calculated tem- 
perature dependence is linear, i.e. the ground state term 
dominates the temperature dependence of F™ a . The ab- 
solute magnitude of F™ a is by about three orders of mag- 
nitude smaller than what we find for FJ 10 and we find that 
we can neglect this contribution. 

The more complicated dependence of the Fermi-Dirac 
factor does not allow us to split fp(l — f~/\ into parts 

as in eq. ©, and the low-temperature dependence of 
F™ is determined by the leading terms in a Sommer- 
feld expansion. The T-dependence therefore is given by 
eq. (|10fl . Studying the sign of the single factors in F™ 
we find that this contribution is positive for all temper- 
atures. Starting from a non-zero constant at T = the 
magnitude slowly approaches zero as the temperature in- 
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creases. Our numerical results confirm this analysis, cf. 
Fig. The constant C2 is the non-adiabatic correc- 
tion to the adiabatic ground state energy*i, &q(V). In 
general C2 <C $0 and is negligible. The curvature A 2 
represents the leading electron-phonon correction to the 
low-temperature free energy and specific heat. The lat- 
ter can be written as 



C = —N bs k 2 B {\ + \)T, (19) 

where the electronic density of states given by band- 
structure alone, Nf, s , comes from 7i e i and is corrected by 
the electron-phonon interaction with the factor A. This 
factor is measurable in experiments and can get quite 
large for certain metals. In the case of lead for instance, 
A w 1.2. We calculated the constant C2 and the cur- 
vature A2 for all four metals we studied and our results 
arc listed in Table ITTT1 In addition to our results we also 
included results from experiment and from previous the- 
oretical studies. Our results agree fairly well with the 
other numbers. It should be pointed out that while F™ 
is important for low temperatures, it is negligible at high- 
T. Comparing the entropy from the non-adiabatic con- 
tribution with the electronic contribution for instance, 
the ratio / S e i while important at low temperatures, 
becomes negligible for T/Te > 1. 

Our results in Tablc llTTI indicate variations in the differ- 
ent theoretical calculations. But regardless of the details 
of the theoretical description, they show the ability of 
the theory to calculate the experimental results. Since we 
are working within the Einstein approximation, they also 
provide confirmation of our previous assumption that the 
Einstein model will give us the correct temperature de- 
pendence and accurate results. 



taken by Wallace for instance^ (Section 25). We used a 
linked cluster expansion, given by 



(Al) 

where the self-energy and Green's function implicitly de- 
pend on the coupling constant™, 77. Since we are using a 
plane wave basis and a free electron dispersion, we need 
to include the zeroth order term in Ti. ep . The electron 
self energy, E(fc, ifc n ), includes diagrams of the form 



U(Q) U(Q) 
► • • ► 

(A2) 

in n th order. We therefore do not consider these terms 
in TL ep . We do need to correct our plane wave basis and 
include terms of the form 




U(Q) h(Q) 

(A3) 

This can be achieved in regular perturbation theory by 
correcting the wave function to first order and then cal- 
culating the diagonal matrix elements in first order. Eq. 
(25.17) of Wallace^ corresponds to this kind of diagram. 
The first term in (25.21) therefore is given by 



APPENDIX A: TECHNICAL DETAILS 
REGARDING THE CALCULATION 

We used a diagrammatic technique to calculate the 
electron-phonon contribution eq. (J2J. To leading order 
in the interaction however this approach is completely 
equivalent to a standard perturbation theory approach 




g(k,Q) 



g(k,Q) 



(A4) 
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with Q 7^ the electronic energies 7^ e^ + g- 
To remind the reader we would like to note that as pointed 
out previously, the first non-adiabatic contribution, F" a , 
and consequently its related entropy, 5™ a , is negligible 
compared to the other two contributions. 



